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ABSTRACT 

In this paper we examine various issues closely related to the ongoing discussion on 
the nature of the Blandford-Znajek mechanism of extraction of rotational energy of 
black holes. In particular, we show that switch-on and switch-off shocks are allowed 
by the shock equations of relativistic MHD and have similar properties to their New- 
tonian counterparts. Just like in Newtonian MHD they are limits of fast and slow 
shock solutions and as such they may be classified as weakly evolutionary shocks. The 
analysis of Punsly's MHD waveguide problem shows that its solution cannot have the 
form of a traveling step wave and that both fast and Alfven waves are essential for 
generating the flow in the guide. Causality considerations are used to argue that the 
Blandford-Znajek perturbative solution is in conflict with the membrane paradigm. 
An alternative interpretation is presented according to which the role of an effective 
unipolar inductor in the Blandford-Znajek mechanism is played by the ergosphcric 
region of a rotating black hole. Various implications of this are discussed. 
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In coordinate systems singular at the black hole horizon, 
like the popular Boyer-Lindquist coordinates, the horizon 
is inevitably turns into a rather peculiar boundary of spa- 
cial domain. According to the widely accepted "Membrane 
paradigm" the horizon, or rather somewhat less stringently 
defined "stretched horizon" which is placed somewhere just 
above the real horizon, may be identified with a rotating 
conducting sphere (e.g. Blandford 1979, Thorne et al. 1986). 
This makes magnetized black holes look analogous to mag- 



netized neutron stars. For many years, beginning with ( Pun- 
sly & Coroniti 1990a), Brian Punsly have been criticizing 



this view and also the perturbative steady-state electro- 
magnetic wind solution for a force-free magnetosp here o f 
a rotating black hole due to Blandford and Znajek (1977), 



the BZ solution, together with similar MHD models (e.g. 
Phinney 1982,1983) on the basis of causality arguments. In- 
deed, in the case of pulsars there is only an outgoing wind 
which passes first through the Alfven critical surface and 
then through the fast critical surface and, thus, the neu- 
tron star can communicate with the wind by means of both 
fast and Alfven waves (Since the gas pressure is dynami- 
cally insignificant in the tenuous magnetospheres of neutron 
stars and black holes, the slow waves seem to be irrelevant.) 
Black holes, however, must also develop an ingoing wind 
which passes through its own pair of critica l surfaces before 
reaching the horizon ( Takahashi et al. 1990 ) . For the typical 



parameters of astrophysical black holes the inner fast surface 
is likely to be extremely close to the black hole horizon and 
one may argue that the stretched horizon can communicate 
with the outgoing wind by means of fast waves. On the other 
hand, the inner critical Alfven surface may be rather distant 
from the black hole horizon. Thus, even the stretched hori- 
zon cannot communicate with the outgoing wind by means 
of Alfven waves which makes black holes rather different 
from neutron stars. 

Blandford (1979) proposed that the outgoing wind of 
black holes is established by means of fast waves alone 
(This seems to be the only way to reconcile the mem- 
brane paradigm with the BZ solution.) However, Punsly 
(1996,2001) argued that fast waves are completely irrele- 
vant and Alfven waves are solely responsible for creating 
the global system of poloidal electric currents of such winds 
and adjusting the angular velocity of magnetic field lines 
and suggested that the steady-state BZ solution is unsta- 
ble and, hence, nonphysical. In partic ular, Punsly c riticized 
Znajek's horizon boundary condition (Znajek 1977) used to 
determine the wind constants in the BZ electrodynamic so- 
lution and in the MHD analysis of Phinney (1982,1983) as 
a condition imposed in a region causally disconnected from 
the outgoing wind. 

However, Znajek's condition is not really a boundary 
condition as it simply prohibits infinitely strong electromag- 
netic field as measured by a physical observer, e.g. a free 
falling observer (e.g. Phinney 1983). Recently, Beskin and 
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Kuznetsova (2000) stressed this point once more and ar- 



nerowicz ( 1967) it was apparently show n tha t limit shocks do 



gued t hat, though the interpretation of the stretched hori- 



not exist in rela tivistic MHD. In Anile (1989) and (Majorana 



zon as a unipolar inductor is misleading, there is nothing 
wrong with causality of the BZ-like MHD models. This con- 
clusion is strongly supported by the res ults of recent time - 
dependent electrodynamic simulations (Komissarov 2001) 
which indicate asymptotic stability of the BZ solution (Zna- 
jek's boundary condition was not imposed in these simula- 
tions.) 

In this paper we continue the discussion of the causality 
paradox a bit further and attempt to clarify the nature of 
BZ mechanism. At first sight, the study presented here may 
appear rather unfocused but a closer look reveals strong con- 
nections between its sections. In Sec. 2 we study limit shocks 
of relativistic MHD. In Sec. 3 we use these results to anal- 
yse the Punsly's waveguide problem that provides important 
insights into the problem of relativistic MHD and electro- 
dynamic winds. In Sec. 4 we propose a modification of the 
membrane paradigm that does not conflict with causality. 



2 LIMIT SHOCKS OF RELATIVISTIC MHD 

It is well known that shock equations of Newtonian MHD 
allow compressive solutions that have non- vanishing tangen- 
tial component of magnetic field, B t , only on one side of 
the discontin uity. These are called swit ch-on and switch-off 
shocks, e.g. (Jeffrey and Taniuti 1964). Both shocks prop- 
agate with Alfven speed relative to the state with non- 
vanishing tangential component of the magnetic field. Rela- 
tive to the other state, a switch-on shock is super-fast and 
a switch-off shock is sub-slow. It is also known that any 
switch-on solution can be considered as a fast shock in the 
limit Bt — > and any switch-off solution as a similar limit of 
a slow shock solutio n and this is why these s hock are often 
called limit shocks (Jeffrey and Taniuti 1964). 

Strictly speaking, switch-on shocks are not evolution- 
ary simply because it is impossible to ensure that B t is 
exactly zero upstream of the shock. An infinitesimally small 
perturbation of the upstream state resulting in an infinites- 
imally small upstream tangential field will generally turn 
this shock into a fast shock followed by an Alfven discon- 
tinuity. This fast shock will be infinitesimally close to the 
original switch-on shock in all respects except the direction 
of the downstream tangential component of magnetic field. 
Moreover, both the fast and the Alfven waves will have in- 
finitesimally close wave speeds. Only a finite amplitude per- 
turbation can result in a finite split of a switch-on shock. 
In this sense, switch-on shocks are similar to fast shocks 
but rather distinct from intermediate shocks of MHD which 
are genuinely non-evolutionary and split as a result of in- 
teraction with waves of infinitesimally small amplitude (e.g. 
Landau & Lifshitz 1959, Falle and Komissarov 2001). Simi- 
lar arguments apply to switch-off shocks. All these specific 
properties of limit shocks explain why they should be con- 



sidered as physically meaningful solutions and why Jeffrey 



and Taniuti ([1964 ) called them weakly evolutionary. 



Shock solutions of relativistic MHD have been a sub- 
ject of rigorous analysis begi nning with the pioneering pa- 
per by Hoffman and Teller (195C). The results have been 
summarized in two r ather comprehensive mono graphs by 
A.Lichnerowicz Jl967|) and and A.M. Anile (|l98S|) . In Lich- 



and Anile 1987) the analysis of limit shocks is not presented 



and readers are referred to the work by Lichnerowicz. It has 
to b e stre ssed that, as it has been explained in Lichnerow- 
icz (1967), the non-existence of limit shocks in relativistic 



MHD would make this system qualitatively different from 
Newtonian MHD. However, the author of this paper has re- 
cently found such shocks in numerical solutions of certain 
Riemann problems and realized that there must be a flaw in 
the analysis of Lichnerowicz. 

Since the style adopted in these monographs is a bit too 
mathematical we shall f ollow the more traditional approach 
of Hoffman and Teller ( 1950 ). First, we construct the limit 
shock solutions, then we show that there exist evolutionary 
shock solutions in the neighborhood of these soluti ons, a nd 
point out the error in the analysis of Lichnerowicz (1967). 

Let us define more precisely what is meant by the limit 
shocks in relativistic MHD. The first condition is that, in 
the rest frame of the fluid on one side of the shock the mag- 
netic field is parallel to the shock normal. This allows us to 
construct a shock frame that is moving with respect to the 
original fluid frame along the magnetic field. In this frame 
the shock is at rest, the electric field is vanishing and the 
fluid velocity is parallel to the magnetic field on both sides 
of the shoc k. Th is frame was first introduced by Hoffman 
and Teller (1950) to simplify the analysis of oblique MHD 
shocks. By construction, in this frame the magnetic field on 
one side of a limit shock is normal to the shock front. The 
second condition is a non- vanishing tangential component of 
magnetic field on the other side. Thus, the shock can be de- 
scribed as either switching on or switching off the tangential 
component of magnetic field in the Hoffman- Teller frame. 

From Maxwell equations it follows that in the Hoffman- 
Teller frame the electric field vanishes on both sides of the 
shock and, thus, the fluid velocity vector is pa rallel to the 
magnetic field vector (Hoffmann fc Teller 1950): 



u x = sB\ (1) 

where u ! a re the space components of the 4-velocity vector. 
Following ( Hoffmann & Teller 195C ) we introduce Cartesian 



coordinates such that the x-axis is along the shock normal 
and B z = on both sides of the shock. Then the only shock 
equations we have to analyse are 



Continuity equation: 
[D] = 0, 

Energy equation: 
[T tx ] = 0, 

Momentum equations: 
[T xx ] = 0, 

[T xy ] = 0, 

where for any quantity A 

[A] = A 2 -A 1 , 

D — pu 



Tti t i 

— wu u , 



(2) 

(3) 

(4) 
(5) 

(6) 
(7) 
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Figure 1. Switch-on shock solutions. Here / = uf/u ai and k = u ai /u si . Index "1" refers to the state with vanishing B v 
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Figure 2. Switch-off shock solutions. Here / = uf/u ai and k = u ai /u si . Index "1" refers to the state with vanishing B v 



T 13 = wu l u 3 + (p + B 2 /2)g l] - B i B j 



(8) 



Here p is the rest mass density, w is the relativistic enthalpy, 
p is the thermodynamic pressure, u v is the fluid 4-velocity, 
and we are using units such that 47r, the magnetic perme- 
ability and the speed of light do not appear in the equations. 
From (UUO) one has 



[{ws 2 - l)B x B y ] = 0. 



(9) 



Thus if B v — only on one side of the shock then on the 
other side 



s 2 = l/w. 

Combining @ with (Q) one obtains 

(n*) 2 = (w + B 2 )/w, 

(u*) 2 = (B x ) 2 /w, 
and, finally, 

(v x f = c 2 a = (B x ) 2 /(w + B 2 ). 



(10) 

(11) 
(12) 

(13) 



The last equation tells us that relative to the state with 
non-vanishing tangential field the shock propagates with 
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Figure 3. Splitting of a switch-on shock into a fast shock and Alfven wave as the result of interaction with a small perturbation of 
magnetic field introduced upstream of the initial solution. The top panels show the initial switch-on solution, where p is the gas pressure. 
The shock is stationary relative to the computational grid. Its internal structure is due to artificial viscosity. The bottom panels show 
the result of interaction with a small perturbation of magnetic field (B z = 0.05 for x > 0.5) upstream of the initial shock. The fast shock 
is located at x ~ 0.05 and the Alfven wave in —0.3 < x < —0.1. 



Alfven speed, c a , e.g. ( Komissarov 1999 ). Substitution of 
these results into (H-M) leads to 



B +2p = 2T X 



w — 2p ■ 



rptQ 



and 

w \B*J ' 



2T X 



(14) 
(15) 

(16) 



where the shock invariants D, T tx , and T xx may be evalu- 
ated given the parameters of the state with vanishing B y . 



Combined with the equation of state, (|14j|15j) allow us to 
determine the thermodynamic parameters on the other side 
of the shock. In the following we use index "1" for the state 
with vanishing B y and index "2" for the state with non- 
vanishing B v . 



It is easy to verify that if v xi ■ 
vanishes. Consider weak shocks with 



u 2 xi = U ttl (l + a) 



then the shock 



(17) 



where u a = c / ~ c a) may loosely be called the "Alfven 
4- velocity" . Substituting this into (|l4-ll6|) and retaining only 
terms linear in a one obtains 
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B£,+2[p] = a(2t0iO, 



[w] — 2[p] = awi, 



(18) 
(19) 
(20) 



perturbation of the upstream magnetic field on a switch-on 
shock. The initial solution 



From these and the second law of thermodynamics one has 
\p] = awiul 1 , (21) 

i), (22) 



B V2 = 2w ^ a ( U ll 



where u 3 = a/yl — a 2 is the "sound 4- velocity", a is the 
sound speed. From ( p2j ) one can see that weak limit shocks 
exist only if 

a(c ai - ai) > 0. 
This can only be satisfied in the following two cases 



ai then p2 > pi , this is a switch- 
then p2 < pi, this is a switch- 



• Ifc ai > a\ and \vf \ > c, 
on compressive shock, 

• If c ai < ai and \vf\ < c, 
off compressive shock. 

Figure [l] shows the switch-on solutions for the poly- 
tropic equation of state with 7 = 4/3, pi = 1, pi = 1 and 



k = ti ai /us-, = 1.1 , 2, 3, 10. Like in Newtonian MHD ( Jeffrey 



and Taniuti 1964), a relativistic switch-on shock turns into 
a pure gas dynamical shock propagating along the magnetic 
field if the shock speed exceeds a certain critical value. 

Figure ^ shows the switch-off solutions for the poly- 
tropic equation of state with 7 = 4/3, pi = 1, pi = 1 and 
k = u a Ju sl = 0.8,0.7,0.5,0.2. For k = 0.8 the switch-off 
shocks turns into a pure gas dynamical shock. For other val- 
ues of k the shock curve terminates as the gas pressure of 
the upstream state vanishes. 

The evolutionary conditions for shock solutions are of 
very general nature and apply to relativistic shocks in ex- 
actly the same manner as to Newtonian shocks. Thus, evolu- 
tionary compressive relativistic MHD shocks must be either 



super- or sub-Alf venic on both sides of the shock (Falle and 



Komissarov 2001). Let us show that in the neighborhood of 
limit shocks there exist evolutionary shock solutions. From 
(^) one can see that if B y has the same sign on both sides 
of the shock then so does 



H = ws — 1. 

From (|l|) one obtains 

{\ + H)B% 



(23) 



(24) 



w + {1 + li)B 2 ' 
Since f'(p.) > 0, one immediately concludes that 

• If fx > then v%. > (? x on both sides of the shock, 

• If n < then v\ < c 2 . on both sides of the shock, 

and, thus, the shock is evolutionary. For a switch-on shock 
fi2 > 0. Thus, if we introduce an infinitesimally small B\ 
of the same sign as B\ the new shock solution will be a 



fast shock, just like in Newtonian MHD ( Jeffrey and Taniuti 



1964). Similarly, a switch-off shock turns into a slow shock. 

Due to the general nature of evolutionary conditions 
splitting of not strictly evolutionary shocks must proceed in 
the same fashion regardless of relativistic or classical nature 
of governing equations. Figure |^ shows the effect of a small 



Left State: 

B = (3.149,2.749,0), 
E = 0, 

p = 3.454, p = 6.743, 



Right state: 

B = (3.149,0,0), 

E = 0, 

p=l, p= 1. 



describes a stationary switch-on shock. It internal struc- 
ture seen in fig. ^| is entirely due to artificial viscosity. The 
upstream state is then perturbed by introducing tangential 
magnetic field with B z = 0.05 upstream of the shock. The 
eventual interaction splits the switch-on shock into a fast 
shock and Alfven wave in the same manner as it has been 
described in Sec.l. In these simulation we used polytropic 
equation of state with the ratio of specific heats 7 = 4/3. 
These simulation were ca rried out using the upwind numer- 
ical scheme described in ( [Komissarov 1999| ). 

Finally, we find nece ssary to explain the error in the 
analysis of Lichnerowicz (1967) that eventually led to the 



incorrect conclusion on non-existence of the limits shock so- 
lutions (p. 161). On page 151 of the book, equation (47-1), 
a space-like 4-vector h a is decomposed into a sum of two 
4-vectors, one parallel and the other one normal to a unit 
space- like vector n°: 

h = t — j]n , t riu — 0. 

Then it is claimed that t a is always space-like. However, 
this is not true. The reader can easily verify that if h a — 
(0, h, 0, 0) and n° = (a/3, 2, 0, 0) then t a is time-like. 



3 PUNSLY'S WAVEGUIDE PROBLEM 

In order to demonstrate the exceptional role of Alfven waves, 
Punsly proposed to consider a much simpler problem involv- 
ing a Faraday wheel connected to a cylindrical waveguide 
uniformly filled with cold tenuous plasma and a strong axial 
magnetic field (Sec. 2. 9. 4 in Punsly 2001). For this problem 
he claimed to have constructed a step Alfven wave solu- 
tion of MHD equations corresponding to an instantaneous 
spinning up of the disc. Downstream of the wave front the 
magnetic field has only axial and azimuthal components, 
the electric field is radial and the electric current is axial 
(see figure ^J) . The return current flows over the surface of 
the guide and the global current closure is ensured by the 
displacement current of the leading front. Such simple step 
wave solution, however, is impossible for the following simple 
reason. 

Since in the unperturbed state the Alfven speed in the 
direction of magnetic field is uniform across the waveguide, 
an Alfven front launched from the surface of the disc would 
indeed stay normal to the guide axis. The problem is that 
such wave would have no effect on the state of plasma, one 
could call it a ghost wave. Indeed, a frame propagating 
with Alfven speed in the axial direction is the Hoffmann- 
Teller frame, as in this frame E — upstream. However, 
in the Hoffmann- Teller frame the amplitude of the tangen- 
tial componen t of magnetic field remains unchanged by an 
Alfven shock (Komissarov 1997) just like it does in Newto- 



nian MHD. Since in this case the tangential component is 
zero upstream it must be vanishing downstream as well. It 
is easy to verify that the Lorentz transformation to the orig- 
inal frame of the waveguide preserves this result. Thus, the 
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Figure 4. The electromagnetic properties of the flow pattern in a cylindrical waveguide attached to a rotating conducting disc as 
described in Punsly (2001). 



Alfven shock has to be ruled out as well as the claim that 
the solution to this problem involves only Alfven waves. 

The tangential component of magnetic can be generated 
by a switch-on fast compression wave including a switch- 
on shock. However, a switch-on shock cannot not remain 
plane given the properties of the "solution" downstream of 
the discontinuity. Indeed, for a constant shock speed rela- 
tive to the unperturbed state the tangential component of 
the downstream magnetic field would have the same am- 
plitude everywhere (see eqs p^jl^ ) whereas in the guide it 
must vanish along the symmetry axis. Thus, the solution to 
the Punsly's problem in the framework of relativistic MHD 
cannot have the simple form of a step wave altogether. 

Although we have not been able to find an analytical 
solution to the waveguide problem, the following basic ar- 
guments suggest that both fast and Alfven waves are impor- 
tant in this problem. When the disc rotation is switched on, 
the induced velocity shear generates the azimuthal magnetic 
field B,/, k r at the disc surface. Since the force balance is 
broken a switch-on MHD shock is driven into the guide. Ini- 
tially this shock is plane but since the shock speed depends 
on r it eventually becomes curved. In the local Hoffman- 
Teller frame of such curved discontinuity the upstream tan- 
gential component of magnetic field will generally have a 
direction which is different from the one of the downstream 
tangential component and, thus, the discontinuity will no 
longer satisfy the MHD shock equations. As explained in 
Sec. 2 it will split mainly into a fast shock followed by an 
Alfven wave. Thus, the disc ultimately emits both fast and 
Alfven waves. Exactly the same evolution would be observed 
in the case of Newtonian MHD. The role of fast waves is to 
create and amplify the non-axial field whereas the Alfven 
waves ensure that it becomes azimuthal. 

In the limit of force-free dege nerate electrodynamics 
( Komissarov 2002 ; Blandford 2002 ) the wavespeed of the 
fast wave tends to the speed of light and this makes a step 
wave solution possible. Indeed, it is easy to verify that equa- 
tions of degenerate electrodynamics allow the following trav- 



eling wave solution: 

B r =0, B 4 ,(t,z) = -Sl(t-z)rB , B z = B , 
E r {t,z)=B 4> {t,z), E^=0, E z = 0, (25) 

where the components of vectors are given in the orthonor- 
mal basis of cylindrical coordinates and Sl(t) is the an- 
gular velocity of the disc. One could be tempted to in- 
terpret the discontinuity of such step-wave solution as a 
limit of a switch-on shock and call it a fast wave. How- 
ever, the eigensystem of degenerate electrodynamics degen- 
erates in the axial dire ction (in general, in the direction of 
E x B ± BVB 2 - E 2 ( |Komissarov 2002[ )), that is both the 
Alfven and the fast linear waves propagate with the same 
speed, which is the speed of light, and are no longer dis- 
tinct. Moreover, both modes are linearly degenerate and, 
thus, their shock solutions have the same properties as the 
small amplitude waves. This makes it quite impossible to tell 
whether this solution represents a fast or an Alfven wave of 
degenerate electrodynamics. In fact, it should be regarded 
as a mixture of both modes. Indeed, in the waveguide which 
is not perfectly cylindrical this degeneracy will be removed 
and the solution will split into two distinct waves in the way 
similar to splitting of a switch-on shock discussed above. 



4 MEMBRANE PARADIGM AND THE 
BLANDFORD-ZNAJEK SOLUTION 

The results of the previous section show that in Punsly's 
MHD waveguide problem both fast and Alfven waves are 
important. This has to be true in general including the 
case of black holes magnetospheres. In such magnetospheres 
Alfven waves propagate only along the poloidal field lines 
(Appendix A). Thus, fast waves must play a special role in 
establishing the cross-field balance whereas Alfven waves are 
particularly important for establishing the wind constants 
(They are constant along the poloidal field lines.) 

Since neither of the waves can be ignor ed the causal- 
ity paradox posed by Punsly and Coroniti ( 1990a ) has to 
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be taken seriously. The resolution of the paradox proposed 
by Punsly and Coroniti involves rejection of both 1) the 
membrane paradigm or rather its part that identifies the 
stretched horizon with the rotating conducting surface of 
a unipolar inductor and 2) the Blandford-Znajek solution. 
They eventually tried to construct completely different mod- 
els of black hole magnetospheres involving rotating dense 
shells and discs of a ccreting mat ter as material analogues 
of the Faraday disc (Punsly 200 f ). However there seems to 



exist a different resolution of this paradox which we discuss 
here. 

In fact, the causality arguments directly hit the mem- 
brane paradigm and it cannot be saved. The analogy be- 
tween the horizon and a rotating conducting sphere is not 
at all that complete as it is believed. It does allow a sim- 
plified presentation for an audience unfamiliar with general 
relativity but it does not really provide deep insights into 
the physics of black hole magnetospheres. 

Things are different when it comes to the BZ solution. 
All what the causality paradox tells us is that this solution is 
inconsistent with the membrane paradigm! Historically, the 
BZ solution was used to build the paradigm and now it is 
widely considered that both are inseparable, but in fact the 
BZ solution clearly indicates the limitations of the paradigm. 
Clearly, there is no material analog of the Faraday disc in 
the BZ magnetosphere. If the horizon does not play its role 
then w hat is forcing the rotation of magnetic field lines? The 



answer |san only lay in the properties of space-time outside of 



the blac k hole horizon as it has been suggested in ( Blandford 
2002|)7J ince another en ergy extractio n mechanism, namely 
the Penrose mechanism ( ^enrose 1969 ), operates only within 
the ergosphere it seems reasonable to consider the possibility 
that it is the ergospheric region of space-time that serves 
as a driving "force" of the BZ mechanism. In fact, this is 
consistent with the causality arguments. 

The super- Alfvenic region of the ingoing wind cannot 
and does not have to communicate with the outgoing wind 
by means of Alfven waves just like the super- Alfvenic region 
of the outgoing wind cannot communicate with the ingoing 
wind in such a way. However, the driving source responsi- 
ble for both the outgoing and the ingoing winds must be 
able to communicate with the winds by means of both fast 
and Alfven waves and, thus, it must be located between the 
Alfven surfaces and, thus, well outside of the horizon (Sim- 
ilar conclusion was reached by Punsly (2001) and Beskin & 



Kuznetsova (200C).) The position of these critical surfaces, 



which merge with the light surfaces in the limit of degenerate 
electrodynamics (Appendix A) , is not fixed as it depends on 
the angular velocity of magnetic field lines, O/, which de- 
pends on many factors including the interaction with the 
surrounding plasma (effective load of the black hole elec- 
tric circuit). However, for all values of fi/ consistent with 
extraction of energy of a black hole in the BZ solution the 
inner Alfven surface is located inside the ergosphere (Ap- 
pendix A). The only exception is the polar direction where 
the light surface, the horizon, and the ergosphere coincide, 
but the Poynting flux density vanishes along the symme- 
try axis and, thus, there is no outgoing wind as well. Thus, 
there always exists an outer region of the ergosphere which is 
causally connected to the outgoing wind of the BZ solution. 

In order to verify this conjecture one could study the 
dynamical behaviour of magnetic field lines that do not pen- 



etrate the horizon. The lines that enter the ergosphere are 
expected to be forced into rotation whereas those that do 
not should remain nonrotating. Such study is under way. 

Naturally, the driving force has to be the same both in 
electrodynamic and MHD models. In MHD approximation 
the wavespeeds of both fast and Alfven waves are smaller 
then the speed of light and therefore the corresponding in- 
ner critical surfaces lay strictly outside of the horizon every- 
where. Since the horizon and the ergosphere always coincide 
in the polar direction, both inner critical surfaces are sit- 
uated outside of the ergosphere in the polar region. Thus, 
there exists a polar flux tube which is causally disconnected 
from the ergosphere. Within such a tube there can only be 
possible an accretion. Such conclusion is not entirely unex- 
pected as a particle located on the symmetry axis is subject 
only to gravitational attraction. 

Another related issue is whether the approximation 
of degenerate electrodynamics (or magnetically dominated 
ideal MHD) brakes down somewhere near the horizon. It 
is impossible to get full answer to this question within the 
framework of degenerate electrodynamics. For example, if 
the condition E ■ B = is satisfied by the initial solu- 
tion it will be preserved during the evolution though there 
may not be enough charged particles to ensure this con- 
dition. However, the preservation of the other condition 
B 2 — E 2 > 0, which is required for the hyp erbolicity of 



degen erate electrodynamics, is not guaranteed (Komissarov 
2002; ). As B 2 -> E 2 / the drift velocity of plasma tends to 
the speed of light indicating that particle inertia may need 
to be taken into account. However, the BZ solution satisfies 
this condition all the way to the horizon and the inertial 
effects are unlikely to be important along the magnetic field 
lines threading the horizon. This is not so obvious in the 
case of magnetic field lines threading the equatorial plane of 
the ergosphere where the approximation of degenerate elec- 
trodynamics or ideal MHD has to break down in order to 
ensure the current closure condition. Punsly and Coroniti 
argue that particle inertia becomes important in 



( |!990b| ) 



the equatorial region as plasma is forced into rotation with 
almost speed of light relative to the zero angular velocity 
observers. However, they have not taken into account the 
Compton drag which may well lead to a much lower Lorentz 
factor of plasma rotation. This problem requires further in- 
vestigation. 



5 CONCLUSIONS 

(i) Switch-on and switch-off shocks are allowed by the 
shock equations of relativistic MHD and have similar prop- 
erties to their Newtonian counterparts. Just like in Newto- 
nian MHD they are limits of fast and slow shock solutions 
and as such they may be classified as weakly evolutionary 
shocks. 



(ii) Contrary to what is claimed in (Punsly 2001), the 



solution to Punsly's MHD waveguide problem cannot have 
the form of a step-like traveling wave and the guide flow 
cannot be established by means of Alfven waves alone. Even 
in the limit of degenerate electrodynamics where a step- wave 
solution exists it involves a mixture both fast and Alfven 
waves. This suggests that both waves are important in the 
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problems of magnetically driven MHD and electrodynamic 
winds. 

(iii) Blandford-Znajek solution contradicts to the mem- 
brane paradigm as the stretched horizon cannot play the 
role of a unipolar inductor. Causality arguments suggest 
that, just like in the case of the Penrose mechanism, the 
driving "force" of the Blandford-Znajek mechanism is the 
ergospheric region of space-time. 
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APPENDIX A: 

The wavespeed of Alfven waves of degenerate electrodynam- 
ics is given by 

(i± = /i ± • n (Al) 
where n is a unit vector normal to the wave front and 



= i f E x B ± Bt/B 2 - E- 



(A2) 



( Komissarov 2002 ) . On a cr itical surface with normal n ei- 

vanishes 



ther /x+ or /i_ given by 
fj, ± ■ n — 



Al 



(A3) 



Consider the Kerr metric in the Boyer-Lindquist coordinates 
{t,<f),r,8}. For a steady-state force-free axisymmetric mag- 
netosphere n has only poloidal component and 



Ftk — —SlfF, 



f r <t>k, 



Ft* = 0, 



(A4) 



where F is the electromagnetic field tensor, Q f is the angular 
velocity of m agnetic field lines, and k — r,9 ( Blandford & 



Znajek 1977). Then in the orthonormal basis {e?, , eg} of 



a local fiducial observer, FIDO, one has 

Ej, = 0, Er = KB^, E § = -KB*, 
where 



K 



(A5) 



(A6) 



a is the "lapse function", and Q,f = —g*t/g<t><j> is the FIDO's 
angular velocity. From A2 , A^ one obtains the poloidal com- 
ponent of n ± : 



Mp ± = {~ KB $ ± VB 2 -K*B$) , (A7) 

whe re B p is the poloidal magnetic field. Thus, the condition 
A3 is satisfied if either 

(B p ■ n) = 0, 

or 

M P± = o. 



(A8) 



(A9) 



A8 simply states that Alfven waves propagate only along the 
poloidal magnetic field lines whereas At is the sought crit- 
icality condition. For B p 7^ this condition can be written 
in terms of Qf and components of the metric tensor as 

f(Sl f ,r, 9) = g^fl} + 2g t *n f + g u = 0, (A10) 



which is the wel l known equation of a light surface (Faka- 
hashi et al. 199fj|) 



Without any loss of generality we may assume that B p 
is outgoing and consider only the northern hemisphere. In 
this case the condition 







(AH) 



corresponds to the inner critical surface whereas 

V p _ = (A12) 



corresponds to the outer critical surface. From (Blandford 



k, Znajek 1977 ) we find that in the case of outgoing energy 
flow 

^L = -B T {l f >0, (A13) 

where £ is the energy flux within a flux tube of magnetic 
flux "I> and 

B T = V=gF rB = a^-B^. (A14) 
Thus, the energy is extracted from the black hole only if 



Btflj < 0. 



(A15) 
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Given this condition, A7 and |A11| show that at the inner 
critical surface 

fi/(fi/ - Of) < 0. (A16) 

For a black hole with positive angular velocity this means 

0<Cl f <n F . (A17) 

Thus, at the inner critical surface the magnetic field lines 
rotate slower than local FIDOs. Similarly one shows that 
at the outer critical surface the field lines rotate faster then 
local FIDOs 

0<Of<%. (A18) 

On the surface of the ergosphere 

f(Qf,r,e)=g^n f (n f -2n F ). (A19) 

From this one can see that in the limit 0/ — > the inner 
light surface coincides with the ergosphere, / being positive 
inside and negative outside. Since 

Jf-=2 flw (n / -nj.)<Q for Q f = 0, 0^0 (A20) 

the inner critical surface moves inside the ergosphere as Qf 
incre ases and must rem ain in side for all values satisfying 
A17 (The third factor in Alt vanishes only when the outer 
critical surface moves inside the ergosphere.) 



